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Summary. In some cases, solutions to nonlinear PDEs happen to be asymp- 
totically (for large x and/or t) invariant under a group G which is not a sym- 
metry of the equation. After recalling the geometrical meaning of symmetries 
of differential equations - and solution-preserving maps - we provide a precise 
definition of asymptotic symmetries of PDEs; we deal in particular, for ease 
of discussion and physical relevance, with scaling and translation symmetries of 
scalar equations. We apply the general discussion to a class of "Richardson-like" 
anomalous diffusion and reaction-diffusion equations, whose solution are known 
by numerical experiments to be asymptotically scale invariant; we obtain an an- 
alytical explanation of the numerically observed asymptotic scaling properties. 
We also apply our method to a different class of anomalous diffusion equations, 
relevant in optical lattices. The methods developed here can be applied to more 
general equations, as clear by their geometrical construction. 

Introduction 



Symmetry methods for the study of differential equations were first introduced 
by S. Lie; they are by now a widely known and very effective set of tools to tackle 
nonlinear equations, both in the sense of theoretical and geometrical study of 
their properties and for obtaining exact solutions. Among the fastly growing 
literature on these, we will here just refer to [20, 27, 37, 42]; see also the shorter 
introductions given in [11, 25, 45]. 

Once an algebra Go of vector fields generating a Lie group Go of symme- 
tries for a differential equation (or system of differential equations) A has been 
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identified, there exist well denned methods to obtain Go-invariant solutions to 
A via a reduced equation A , sometimes denoted as A = A/Go- These are 
also called partially invariant solutions, as in general Go is only a part of the 
full symmetry algebra of A. This method is exposed in detail in most books on 
symmetry methods for differential equations, see e.g. [37] or [45], and here we 
will just briefly recall it. 

Other methods, generalizing or extending this approach, have been devel- 
oped by several researchers [39]; for recent and quick reviews see e.g. [10, 24]. 
Among these we mention in particular the method of conditional symmetries 
and conditionally invariant solutions [13, 18, 28, 41]; and an extension of this, 
the method of partial symmetries 1 [12] (related in turn to "weak symmetries" 
[38]), which will be of use later on in our discussion. 

No matter how successful the symmetry approach, it has long been known 
that there are equations which exhibit asymptotically stable solutions u*(x,t) 
which are invariant under an algebra G* which is not a symmetry of the equa- 
tion; moreover, solutions which are asymptotic to u*(x,t) are in general not 
invariant under G*, so that G* is only an asymptotic symmetry of these. 2 

Several authors have attempted a description of asymptotic invariance (mainly, 
scaling and/or translational) properties of the attracting asymptotic solution 
u*(x,t) mentioned above; it has been generally agreed that such a description 
should be based on (a suitable version of) renormalization group concepts. We 
recall here the books by Barenblatt [4, 5] and by Collet and Eckmann [14], to- 
gether with the works of Bricmont and Kupiaincn [7] and of Goldcnfcld et al. 
[23]. We also quote the rather abstract approach proposed some time ago by 
one of us [20], which will be developed here, focusing on the physically most 
relevant case of asymptotic scaling symmetries, into a concrete theory provid- 
ing an explanation to observed asymptotic symmetry properties in the class 
of equations mentioned above, and possibly prediction of similar properties in 
other systems. 3 

As mentioned above, in this note we will consider scalar PDEs (and focus in 
particular on those of reaction-diffusion type associated to anomalous diffusion; 
see e.g. [9] for a recent review, focusing on aspects of interest here), and accord- 
ingly we will adopt suitably simplified notation and definitions; note however 
that the methods developed here do generalize to the vector case. 

1 The reader should be warned that the terms "partial symmetry" or "partial invariance" 
are also used with a different meaning in the literature, see e.g. [39, 40] . 

2 Among equations exhibiting such a behaviour are a wide class of diffusion and reaction- 
diffusion equations related to anomalous diffusion (studied numerically in [31], to which we 
refer for a discussion of their physical significance), which we will study in this note as a 
nontrivial example of application of our method; these will be called "Richardson-like" as 
they include the Richardson equation for developed turbulence. Another class of equations 
with such a behaviour - also studied below - describes the marginal Wigner distribution 
in certain optical lattices [1, 29], which have recently attracted some attention due to the 
unexpected appearance of Tsallis distribution as stationary states. 

3 We stress that we focus on asymptotics for large t, and in the cases of interest here this 
will correspond to large \x\ as well. Asymptotic - or "approximate" - symmetries for small 
|rr| have been considered by a number of authors, see e.g. [3, 11, 19, 22]. 
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The paper is organized as follows. In sect.l we introduce in some length the 
basic concepts of symmetry for differential equations; and in sect. 2 we recall the 
notions of conditional and partial symmetry. We pass then to new material. In 
sect. 3 we introduce our definition of asymptotic symmetries (discussing directly 
the case of second order scalar evolution equations), and in sect. 4 we extend 
this to conditional and partial asymptotic symmetries. In sect. 5 we discuss how 
these can be used to investigate the asymptotic behaviour of an evolution PDE. 

The rest of the paper is mainly devoted to a detailed investigation of a model 
class of diffusion and reaction-diffusion (RD) equations related to "Richardson- 
like" anomalous diffusion, and their asymptotic behaviour, making use of the 
method proposed here. In sect. 6 we deal with the anomalous diffusion case, 
which turns out to be rather simple. In sect. 7 we discuss the standard FKPP 
equation from the point of view proposed in this paper, discussing in particular 
the asymptotic symmetries of the equation and of its asymptotic solutions. In 
sect. 8 we consider Richardson-like FKPP equations, recall the outcome of nu- 
merical experiments on these, and determine their asymptotic scaling-invariant 
solutions; this provides an explanation of experimentally observed scaling prop- 
erties. In sect. 9 we consider again Richardson-like FKPP equations, but look 
at more general asymptotic symmetries, focusing on large t behaviour near a 
propagating front. It turns out this analysis provides again a direct explanation 
of the asymptotic solutions observed in numerical experiments. 

In section (10) we consider a different anomalous diffusion equation, i.e. 
the one describing the marginal Wigner distribution in momentum space of an 
atom in an optical lattice [29]; attention has been recently called to this not 
only for its relevance in applications, but also because - somewhat surprisingly 
- its asymptotic, time- invariant, solutions are in a certain regime scale-free 
and correspond to a Tsallis statistics [43] . We show that this behaviour can be 
explained on the basis of the asymptotic symmetry properties (under generalized 
scalings) of the equation. 

The final sect. 11 is devoted to summarizing and discussing our findings, and 
presenting some final remarks. 

Acknowledgements. We thank D.Lcvi for useful discussions. The work of 
GG was supported in part by GNFM—INdAM under the program "Simmetrie e 
tecniche di riduzione per equazioni differenziali" ; the work of RM was supported 
by I.N.F.M. (Istituto Nazionale di Fisica della Materia). 

1 Symmetries of differential equations 

In this section we recall the main concepts and definitions for symmetries of 
differential equations (see e.g. [20, 37, 42, 45] for general treatments of these). 
As in this note we will only consider scalar PDEs, we will just specialize our 
formulas to this case. 
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1.1 General notions and properties 

We consider an equation of order n for u = u(x) G R, with x G R m . We 
denote by B = R m the space of independent variables, by U = R the space of 
dependent variables, and by M = B x U the total space of independent and 
dependent variables. 

Note that M has a natural structure of fiber bundle over B, i.e. (M 7 ir,B) 
is a bundle with fiber n~ 1 b — U. A function u = f(x) is then represented by a 
section 7/ of this bundle. 

A vector field in M, written in coordinates as 

d d 

x = eM d^ + v{x > u) d^> (L1) 

will generate a one-parameter group of transformations in M; at the infinitesimal 
level these are described by x l — > x 1 = x' L + s^ l (x, u), u — > u = u + sip(x, u). 

Thus, for e sufficiently small, the function u = f(x,t) is mapped to a new 
function u = f(x,t); one obtains with standard computations that this new 
function is given by 



f(x) = f{x) + e 



• ( ^ ) ' 



+ o(e) (1.2) 

M=/(x) 



where di := d/dx l ; functions depending on u should be computed in u = f(x). 

We introduce some standard notations: first of all, we denote by Dj the total 
derivative with respect to x l , i.e. Di := [d/dx l ) + Ui(d/du) + Uij(dfduj) + .... 
We will use multiindiccs J = (ji, j m ); the order of J is |J| = ji + ... + j m . 
With these, we write Dj := (Dj)- 7 ' 1 ...(D m y m , and uj = Dju; we also write 
uj.i '■— DiUj. Finally, the space of (x, u) and x-derivatives of u up to order k is 
said to be the jet space of order k for M, and denoted by J^M; it is convenient 
to also use the formal limit k — > 00, in which case we write J^M. 

By differentiating (1.2), we get formulas for the transformation of the partial 
derivatives of any order of u with respect to x and t under the action of X. 

The action of the vector field X = £ % di + (pd u on the jet space J^M is 
described by the prolongation of A, i.e. by the vector field 

A* :^ A + (1-3) 

where the coefficients are given by the prolongation formula 

^ = D 3 (^-u a k e)+u a j. k e ; (1.4) 

this reads also, in recursive form and with := <p a , 

= D k *3-u 3ti D k e . (1.5) 

The differential equation A of order n will be written as $(x, u, u^) = 
for some function $ : J^M — > R. We can then apply the vector field to 
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$ and consider the result of this 4 once the constraint $ = 0, i.e. the equation 
A itself, has been taken into account. 

If $ = u t — F(x,t,u,u X7 u xx ), as in our applications below, taking A into 
account consists simply in writing F(x, t, u, u x , u xx ) for u% (this also extends to 
differential consequences: e.g., u xt — D X F). 

We say that X is a symmetry of the equation A given by <f> = if 5 

= . (1.6) 

J $=o 

Condition (1.6) guarantees that X maps solutions to A into (generally, dif- 
ferent) solutions to A; conversely, if X maps solutions into solutions, then (1.6) 
is necessarily satisfied. Actually, the qualitative meaning of symmetries of a 
differential equation is precisely that of transformations mapping any solution 
to A into a solution [21, 37, 42, 45]. 

These concepts can also be stated in an equivalent but more geometric way, 
more convenient for our discussion, as follows. The equation A of order n 
identifies a solution manifold S in J^M; a function u = f(x) is a solution to 
A if and only if the prolongation 7^ of the section 7/ satisfies 7^ C S. The 
vector field X is a symmetry for A if its prolongation X^ to J^M is tangent 
to 5, i.e. 

x (n) . S ^ TS . (1.6') 

We stress that in this note we are mainly interested in two kinds of symme- 
tries: scalings and shifts (or translations). We recall that the shift x % — ► x % +esi, 
u — ► u + ss (with Sj e R) is generated by X = sfli + s a d u , and the scaling 
i x t , u — > X k °u (with ki € R) is generated by X — kixdi + k ud u . 



1.2 Symmetry reduction for PDEs 

Let A be a PDE for u = u(x), x e R p , admitting a vector field X as symmetry; 
in order to determine X-invariant solutions to A one proceeds in the following 
way. (This is standard material, recalled here to fix notation; for more details, 
see e.g. the discussion in chap. 3 of [37].) 

We write A in the form F{x,u^) — 0, with F a smooth scalar function, 
F : J^M — > R, and denote by Y the fc-th prolongation of X. For reaction- 
diffusion equations, p = 2 and k = 2. 

First of all we pass to symmetry-adapted coordinates in M. In practice, 
we have to determine a set of p independent invariants for X in M, which we 
will denote as (y 1 , y p ~ 1 ; v): these will be our X-invariant coordinates, and 
essentially identify the G-orbits, while the remaining coordinate a will be acted 
upon by G. In other words, G-orbits will correspond to fixed value of the (y, v) 
coordinates, with a taking values in a certain subset of the real line. 

4 Note that if we apply to <J>, only terms corresponding to |J| < n in (1.3) will actually 

matter: we can as well consider the (formal) infinite-order prolongation. 

r 'X is a strong symmetry of A if (1.6) holds without the restriction to A (i.e. to $ = 0). 
Any strong symmetry is also, of course, an ordinary symmetry; see [8] for details of the relation 
between ordinary and strong symmetries. 
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The invariants will be given by some functions of the x and u, i.e. 

y l = r ] i (x,u) (i = l,...,p- 1) , v = ((x,u), (1.7) 

with a = a(x, u) as well. Assume (this is always the case for A a scaling vector 
field) that we can invert the above for x and u as functions of (y, v; a): 

x % = x'(y,w;cr) (i = 1, ...,p) , u = 0(y,v;a) . (1.8) 

If now we decide to see the (y; a) as independent variables and the v as the 
dependent one, we can use the chain rule to express ^-derivatives of u in terms 
of the a and y-derivatives of v. Using these, we can finally write A in terms of 
the (a,y;v) coordinates, i.e. in the form $(ct, y, w^ n '). As A is a symmetry of 
A, it follows that the function when subject to the side condition dv/da = 0, 
is independent of a. 

In fact, dv/da = expresses the fact that the solutions v = v(y,a) are 
required to be invariant under X: the equation obtained in this way represents 
the restriction of A to the space of G-invariant functions, and is therefore also 
denoted as A/G. 

Suppose we are able to determine some solution v = f(y) to the reduced 
equation; we can write this in terms of the (a;, it) coordinates as 

C(x,u) = f[ V (x,u)} (1.9) 

which yields implicitly u = f(x), the corresponding G-invariant solution to the 
equation A in the original coordinates. 

Remark 1. Given any vector field X we can consider its characteristic Q := 
<f — Ukl; k and its evolutionary representative Xq := Q(d/du). By (1.4) or (1.5), 
the prolongation of the latter will simply be Xq = Xq + ^2j(DjQ a )(d/duj), 
and standard computations [37] show that X* = Xq + £ k Dk. These formulas 
give a very convenient way of computing prolongations. 

The reduction procedure described above can also be described in a slightly 
different way using evolutionary representatives: if we look for A-invariant 
solution u = f(x) to A, we determine the characteristic Q = ip — u^ 1 of 
the vector field X, and supplement A with the equations DjQ = with 
|J| = 0, ...,k — 1. The equations Q = require that the evolutionary repre- 
sentative Xq = Q(d/du) of X vanish on jf, i.e. that u is A-invariant, and 
all the equations with | J| > are just differential consequences of this. The 
A-invariant solutions to A are in one to one correspondence with the solutions 
to the system Ax := {A; DjQ — 0}. See e.g. [45] for details, and for how this 
approach is used in a more general context. 

Remark 2. The standard method discussed here applies under a nondegeneracy 
(transversality) condition, guaranteeing a certain Jacobian admits an inverse. 
When this is not the case the treatment should go through the approach de- 
veloped by Anderson, Fels and Torre [2]; in case of partial transversality see 
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also [24]. Also, this method is justified only if the (possibly, only local) one- 
parameter group generated by A has regular action in M. Both these conditions 
are satisfied for scaling vector field other than the trivial A = ud u (which is a 
symmetry only for linear equations), and translations other that X — d u (which 
is a symmetry for equations of Hamilton- Jacobi type). 

1.3 Solution-preserving maps 

Consider a PDE A, say of the form 

A ee u t - F(x,t,u,u x ,u xx ) = 0, (1-10) 

with {x,t,u) e M = R 3 , and M a bundle (M, n, B) over B = R 2 . Let us 
now consider another manifold M — R 3 which is also a bundle (M, n, B) over 
B = R 2 , and a C n map F : M — > M mapping (x, t,u) £ M into £, t, w) £ M, 
with of course £ = £(x, t, u), r = r(x, t, u), w = w(x, t, u)). This can be lifted to 
a map once we identify w as the dependent variable in 

M. In this way, (ut,u x ;u t t,u x t,u xx ) are also mapped to expressions involving 

(w T , W(] W tt , W£ T , W£t)- 

The geometrical meaning of this is that the solution manifold S C J^M of 
A is mapped by fW into the solution manifold S C J^M of A. 

If the map is projectable [21], i.e. such that r = r(i), £ = £(ar, t), then (1.10) 
is mapped into an equation of the same type: 

A EE W T — G(£,T,W,W£,W££) = . (1-11) 

Recall now that the function u = f(x,t) corresponds with the section 7/ = 
{(x, t, f(x,t))} in (M,n,B). If F is projectable, we are guaranteed that the 
manifold 7/ C M is mapped into a manifold 7 g C M which is a section for 
(M,9,B). In other words, F(7/) = 7 9 = {(£,T,ff(£,T))} C M. This obviously 
extends to prolongations, i.e. : 7^ — > jg n \ This guarantees, in particular, 

that solutions u = f(x, t) to A are mapped into solutions w = r) to A. We 
say therefore that F is a solution preserving map [44]. 

If F is invertible (with a C™ inverse) , we can repeat these considerations for 
^> = F~ x . In this case, therefore, the equations A and A are equivalent, in that 
there is a C n isomorphism between solutions to A and solutions to A. 

2 Conditional and partial symmetries 

As mentioned above, A is a symmetry if it maps any solution to a (generally, 
different) solution. However, there will be vector fields which map some solutions 
into solutions, and some other solutions into functions which are not a solution 
to A. In this case we speak of partial symmetries. 

That is, A is a partial symmetry for A if there is a nonempty set Sx of 
solutions to A which is mapped to itself by A. Note that Sx could be made 
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of a single solution, and more generally that there could be solutions which are 
left invariant by X. In this case, we say that X is a conditional symmetry for A. 
Obviously any symmetry is also a conditional symmetry, and any conditional 
symmetry is also a partial symmetry. 6 

It is clear from (1.2) that u = f(x) is invariant under X if and only if 

A x := <p[x, f(x)] - u£[x, f(x)] = ; (2.1) 

thus A-invariant solutions to A are solutions to the system A of differential 
equations made of A and Ax- In other words, X is a conditional symmetry of 
A if and only if it is a symmetry of A. 

Partial symmetries can be seen in a similar way. As discussed in [12], a 
function / in the globally invariant set of solutions to A : F = 0, / e Sx in the 
notation used above, will be a solution to a system 

F (0) = o 

^ (1) =° (2.2) 
I = 

where = F, and F( fe+1 ) := X*[F^]. The integer p, i.e. the order of the 
system, is determined as the lowest order such that F^ vanishes identically on 
solutions to the previous equations. Note also that each equation F^ — can, 
and should, be simplified by taking into account the previous equations. 

Needless to say, if A and A are related by a solution-preserving map <&, this 
entails a relation between their respective conditional and/or partial symmetries. 

For a recent discussion on partial symmetries and related notions, see [10]; 
for examples, see [12]. 



3 Asymptotic scaling symmetries for evolution 
PDEs of second order 

As suggested by their names, asymptotic symmetries of a differential equation 
A are vector fields X which, albeit in general not symmetries of A, satisfy 
X^ : Sa — * TSa asymptotically (see below for the precise sense of this). Any 
exact symmetry is also a (trivial) asymptotic symmetry. 

Asymptotic symmetries were introduced and discussed in quite a general 
framework in [20], see also [11]. In this section we introduce them in a simpli- 
fied scheme, i.e. limiting our considerations to the class of equations of interest 
here (evolution PDEs, in particular of reaction-diffusion type) and to the trans- 
formations we are interested in (scalings and shifts). 

6 Partial symmetries (as defined here) were introduced and characterized in [12]. The origin 
of conditional symmetries could maybe be traced back to Felix Klein [20]; for their modern 
theory, sec [13, 18, 28] (here we adopt the point of view of [28], see also [45]). For the relation 
of conditional symmetries to the "direct method" of Bluman and Cole [6], see [41]. 



8 



3.1 Asymptotic symmetries for equations of RD type 



Thus, we discuss scalar equations A for a real function of two variables, u{x, t). 
In the notation of sect.l this means B — R 2 , M = R 3 . We write the equation 
A in the form u t — F(x, t, u, u x , u xx ) = 0. 

We denote by T the space of maps F : (x,t,u,u x ,u xx ) — > R which are 
polynomial in (u, u x , u xx ) and rational in (x,t). Note that this space, which 
corresponds to the space of equations in the form we are considering, is invariant 
under scaling transformations and under translations. 

Vector fields will be written as X = rdt + £,d x + tpd u , and we mainly restrict 
to scaling ones, i.e. 

X = (at)d t + (bx)d x + {cu)d u . (3.1) 

This generates the one-parameter group of scaling transformations 

t -» X a t , x -» X b x , u — > X c u . (3.2) 

The second prolongation of X, which we denote for ease of writing as Y = 
X^ 2 \ is again a scaling vector field, 

y = X + [(c - a)u t ]d Ut + [(c - b)u x ]d Ul . . 

+ [(c - 2a)u tt ]d Utt + [(c-a- b)u xt ]d Uxt + [(c - 2b)u xx ]d Uxo: , 

and generates the scaling group 

t — > \ a t , x — » A b x , u — » A c u , 

ut -> A c " Q u t , u x -» A c_fc u :E , (3.4) 
Utt -» \ c ~ 2a u tt , u x t -» A c_a_b u a;t , u xx -> \ c - 2b u x x . 

Under the action of X, the function u = f(x,t) is transformed into u = 
f(x,t) with, see (1.2), 

f(x,t) = f{x,t) + e [u — atu t — bxu x ] + o(e) . (3.5) 

Thus X induces a flow X in the space of functions u = f(x, t), or more precisely 
in the space of sections of the bundle (M,n,B). Note that X corresponds to 
the evolutionary representative Xq of X, see sect.l. 

We write the flow issued by fo{x,t) as f\(x,t) — exp[AX], with df\/d\ = 
X[f\]- We say that fo{x,t) is X-invariant if it is a fixed point for the flow of 
X, i.e. ifX[/ o ]=0. 

It is also possible that f is not invariant under X, but the flow f\ is asymp- 
totic to an invariant function /», i.e. 

lim \f x (x,t) - U{x,t)\ = , X[U] = . (3.6) 

When (3.6) is satisfied, we say that / is asymptotically X-invariant under the 
flow of X. 
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3.2 Symmetries in the space T 

Let us now consider an equation A : t, u^) = written as 

A : u t - F (x,t,u,u x ,u xx ) = 0. (3.7) 

Such an equation is thus identified by the function F , and we assume F e T . 

The action of X on the space of such equations is described by the second 
prolongation Y = X^ 2 \ In order to consider the flow induced by Y on the space 
of equations of the form (3.7), i.e. on T, we write 

A A = e Ay A := cr(A) [u t - F\(x, u, u x , u xx )] ; (3.8) 

by construction, this satisfies 



dX 

With A as in (3.7), 



Y(A X ) . (3.9) 



Y[A] = \ c - a u t - 



\ b dF xc dF , c _ 6 dF ^ c _ 2b OF 

X b x— + X c u— + X c b u x — + \ c 2b u xx - 

ox ou ou x ou xx 



(3.10) 



We wish to identify Aa with the corresponding F\, see again (3.7). Hence 
we also write, comparing (3.7) and (3.10), 

f = X a+b - c x^f + A-«^ + \«-\ x ^- + X^ b u xx ^- . (3.11) 
d\ ox Ou Ou x ou xx 

In this way, X induces (via Y) a vector field W in the space F; we rewrite 
(3.8) as dF/dX = W(F). 

As recalled above, X is a symmetry of A if and only if Y : S — > TS (recall 
Y = X^). This condition is now rephrased in terms of T by saying that X is 
a symmetry of Ao given by (3.7) if and only if Fo is a fixed point for the flow 
of W. 



3.3 Asymptotic symmetries in the space T 

Suppose now that X is not a symmetry of Ao - Fo is not a fixed point for the 
flow of W in T - but that the flow issued by F under W satisfies 

lim \F X -F,\ = , W(F,) = 0. (3.12) 

In this case we say that X is an asymptotic symmetry for F , i.e. for the equation 
A . 

Remark 3. The procedure given here to define asymptotic symmetries can 
be reinterpreted in a slightly different way: that is, we combine the action of 
W ~ Y = X^ with a rescaling in the equation, so to keep this in the form 
u t - F = 0. 
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By construction, and by the invertibility of the scaling transformation (3.4) 
for A finite, if u = f(x, t) is a solution to Ao, then u = f\(x, t) will be a solution 
to A\, and conversely if u = f\(x, t) is a solution to A\, then u = f(x,t) is a 
solution to Ao. 

This is represented in the following diagram, where 5a C J^M is the 
solution manifold for Ay. 

u = fo(x,t) w 70 C 

u\ = f\(x,t) w 7A C 

In the limit A — > oo the invertibility of (3.4) fails; we can nevertheless still 
say that solutions to the original equation flow into solution to the asymptotic 
equation, provided both limits /» and A» exist. 

If the limit A* exists, it captures the behaviour of A for large A, i.e. for 
large (or small, depending on the sign of a and b) t and \x\. 



So 



Sx 



A 



(3.13) 



3.4 Solution-preserving maps and asymptotic symmetries 

If there is an equation A whose asymptotic behaviour is well understood, and 
such that there exists a solution-preserving map $ : A — > A, we can of course 
study the asymptotic behaviour of solutions u(x, t) to A by means of the asymp- 
totic behaviour of solutions u>(£,t) to A. This approach is schematized in the 
diagrams below: 

A o — ► uo — ► wo 

\e™ e>* e»* (3.14) 

A x - — ► Aa u ^ — » ^ 

Here W := <J>l n) (VK) is the push-forward [35, 36] of the map $("), i.e. of the 
lift of $ : Mi -> M 2 to a map : j( n )Mi -> J {n ^M 2 , to the tangent map 
$ (n) . t[J(")Mi] -» T[J(")M 2 ]. 

If $ is a solution-preserving map, these diagrams commute; if $ has a smooth 
inverse we can study the flow of W using 

e XW = [$ (n)]-l Q e AW Q $(n ) _ (3 15) 

Note that (3.15) remains valid also in the limit A — > oo, i.e. for the asymptotic 
regime. 

Needless to say, this approach is particularly convenient when the asymptotic 
behaviour of the solutions w(£, r) to Ao - i.e. the behaviour of solutions to A* - 
is simple and/or in some way universal (even better if Ao is a fixed point under 
W). 

Denoting as tu*(£, r) the limit expression for the solutions to Ao, the asymp- 
totic behaviour of the solution u(x, t) to A will be given by 

= K(£,t)] . (3.16) 
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More precisely, considering the sections 7/ and "/ g corresponding to u = 
f(x,t) and w = r) (sec scct.l), and going asymptotically into 7^ and 7* 
respectively, we have 7^ = <i> _1 (7*). 

4 Conditional and partial asymptotic symme- 
tries 

It is also possible to consider asymptotic versions of conditional and partial 
symmetries. These are introduced in accordance with our general scheme. 

That is, we could have the case where a function u = fo(x, t) (more precisely, 
the corresponding section 70 C M) flows under e xx to a fixed point u = f*(x, t) 
(more precisely, a Jf-invariant section 7, C M) albeit Ao does not flow to a 
fixed point 7 . 

In such a situation, the solution manifold S\ does not go to a limit manifold, 
but there is a submanifold of this, with 7 a C C S\, which flows to a 
fixed limit submanifold , with 7, C . 

In this case we say that X is a conditional asymptotic symmetry for A. 
Examples of this phenomenon will be discussed in sect. 8 below, where we deal 
with generalized FKPP equations. 

The same construction, with suitable and rather obvious modifications, ap- 
plies for what concerns partial symmetries: suppose that A = A^ ' does not 
flow to a fixed point under the vector field W induced by X in T. Consider the 
equations 

A« :=y[A(°>] , ... , AW :=Y[A^} 

up to an r if it exists - such that A' r ^ does admit a fixed point A^ under 
the W flow, while the A( fe ) with k < r do not. Then the manifold 

s{ 0) n...ns{ r) := s x 

(with C j( n ^M the solution manifold for A^) flows to a limit submanifold 
5*, and solutions u = fo{x, t) to the system A( fc ) (k = 0, r) flow to functions 
u = f*(x,t) such that the prolongation 7*"^ C J^M of the corresponding 
section jf t = {(x, t, f*(x, t))} lie in 5*. In this case we say that X is a partial 
asymptotic symmetry for A. 

Remark 4. If A and A are related by a solution-preserving map <&, this also 
relates their conditional and partial symmetries. In particular, if $ : A — > A is 
invertible and A admits X as a conditional symmetry, then A admits < &~ 1 (X) 
as a conditional symmetry, with the push- forward of This will be of 

use below. 

7 In the general framework [20] the situation would be more complex, as the invariancc 
vector field X would not in general be the vector field applied to extract the asymptotic 
behaviour. Here, however, we do not need to worry about such a general setting. 
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5 Asymptotic symmetries as a tool to test asymp- 
totic behaviour 



In many physically relevant cases, one has to study nonlinear PDEs which are 
not amenable to an exact treatment, or at least for which such a treatment is 
not known, and for which a numerical study shows an asymptotic behaviour 
which appears to be well described by some kind of invariance. Usually, the 
latter corresponds to a scale invariance (self-similar solutions), or a translation 
invariance (travelling waves), or a combination of both of these. 

The discussion conducted so far can be implemented into a procedure allow- 
ing on the one hand to test if the observed asymptotic behaviour is a charac- 
teristic of the equation (rather than an artifact of the numerical experiments 
conducted on it), and on the other hand to formulate simpler equations extract- 
ing the asymptotic behaviour. 

We will now describe the procedure in operational terms; as a (rather simple, 
but relevant) example to illustrate our procedure we will consider here the heat 
equation, while in later sections we apply the procedure on equations associated 
to anomalous diffusion. 

We will denote by X the vector field describing the observed invariance. For 
the sake of concreteness, we consider a second order equation for u = u(x,t) of 
the form u t — F(x,t,u,u x ,u xx ). We denote, as usual, by M the space R 3 of 
independent and dependent variables, i.e. M — {(x, t;u)}. 

• Step 1. Pass to symmetry-adapted coordinates in M, i.e. coordinates 
(a, y;v) such that X(y) = X(v) = 0; thus in these coordinates X = 
f(a,y,v)(d/da). 8 

• Step 2. Identify v as the new dependent variable, i.e. v = v(a,y). This 
allows to write x and t derivatives of u as a and j/ derivatives of v, hence 
to write the differential equation A(x,t;u^) as A(<r, y; v^). 

• Step 3. Reduce the equation A(cr, y; v^ ) to the space of X-invariant 
functions, i.e. to v(a,y) satisfying v a = 0. For X an exact symmetry, the 
reduced equation Ax will not depend on a at all. For X a conditional 
or partial symmetry, a will still appear parametrically in the reduced 
equation. 

• Step 4. Study the asymptotic behaviour of the solutions to the reduced 
equation Ax for a — > oo. 

• Step 5. Go back to the original variables. 

8 We stress that we are not requiring / = 1; actually when we deal with scaling symmetries 
it is appropriate to require /(<r, y, v) = <r. 
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An elementary example: the heat equation 

Let us illustrate our procedure by applying it on the heat equation u t = u xx . 
Its asymptotic solutions arc of the form 



u(x,t) = —= exp 
V2t 



Ax 2 

and are invariant under the scaling vector field 

X = xd x + 2td t - ud u , 



(5.1) 



(5.2) 



which is also an exact symmetry of the equation. 

The symmetry-adapted coordinates (and the old ones in terms of these) are 



a = t , 
t = a , 



y = x 2 /t , v = xu ; 
x = \f&y -i u ~ v/a . 



Using these, the heat equation reads 

av a = 4yv yy + (2 + y)v y + v 



(5.3) 



(5.4) 



Imposing v a — we obtain an equation with a disappearing completely. Need- 
less to say, the equation obtained in this way has solutions 



v{y) = A^/y exp[-y/4] , 



(5.5) 



which when mapped back to the original coordinates produce the gaussian (5.1). 

Note also that (5.4) necessarily requires that for a — > oo, the function v = 
v(<j,y) satisfies v a = 0: thus, the full asymptotic behaviour of (5.4) is captured 
by (5.5). 



6 Asymptotic symmetry of Richardson-like anoma- 
lous diffusion equations 

In this section we apply our procedure to anomalous diffusion (or Fokker-Planck) 
equations 

u t = L[u] (6.1) 



with L the linear operator 



L\u] := 



dx 



t 1 



x x - a l 2 u. 



(6.2) 



These will be called "Richardson-like" (in the following, RL) since for a 
2/3 and v — 3/2 wc get the Richardson equation 



c 2/3 — 
dx 



x 2 ' z u, 
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describing the evolution of the distance between two particles in developed tur- 
bulent regime. 

Needless to say, equations in the class (6.1), (6.2) are much more general. 
We mention that for a = 2 we have the generalized gaussian process 

Ul t Hxx • 

We have the following general result, confirmed (and actually suggested) by 
numerical experiments on a number of RL equations [31]. 

Theorem 1. The universal asymptotic solution to RL equations (6.1), (6.2) 
for L 1 initial data is given by 

u(x,t) ~ Icxph^A 2 "] . (6.3) 



Proof. In the general case described by (6.1), (6.2), the change of variables 

r = t av , £ = x a ' 2 , w = tS 2 - a ^^u (6.4) 

maps the equation (6.1), (6.2) into the heat equation 

w T = vj££ . (6.5) 

Using the inverse change of coordinates we have that the universal asymptotic 
solution 

w(£,t) ~ t~ 1/2 e-« 2/r (6.6) 

of the heat equation is mapped back into the function given in the statement. 
This represents therefore the universal asymptotic solution of the equation (6.1), 
(6.2). 

Remark 5. It is appropriate to stress that this result applies to a quite large 
class of stochastic processes. This includes processes with non-independent but 
finite variance increments (coloured noise), and processes with independent in- 
crements but infinite variance (Levy flights); see e.g. [30, 31] again for examples 
and detailed numerical studies. 

The relevant point is that, in order to apply our approach, the consid- 
ered stochastic processes should present an asymptotic probability distribution 
(starting with sufficiently localized initial distribution) of the form 

P(x,t) ~ pt~ u Q(x/t u ) (|a;| — > oo , t — > oo) , (6.7) 

with p a normalization constant. E.g., for coloured noise Q(z) = e~ z " . 

Remark 6. We also note that one can use this approach for stochastic processes 
which are not described in terms of a diffusion process, but which show in 
numerical experiments a limiting distribution of the form (6.7). This is the case 
e.g. for the random walk on a comb lattice (in this case v = 1/4, and Q(z) = 
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exp[— z a ], a = 4/3); or for Levy flights (here a < 3, v = (a — > 1/2 and 
Q(z) = z~ a ), see [31] for details (the equation studied there in relation to Levy 
flights is a phenomenological one, corresponding to the propagator). In these 
cases we can use our argument to provide an effective equation for the limiting 
behaviour of the probability distribution; this will provide, by construction, the 
correct behaviour in the large \x\ and t limit. 

7 The FKPP equation: asymptotic solution, and 
symmetries 

We now apply our approach to the Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) 
equation [17, 26] reads 

u t = Du xx + eu(l - u) , 

with e and D real positive constants. Here D represents a diffusion coefficients, 
while the parameter e is the inverse of the timescale for a logistic growth. We 
are interested in solutions such that u(x, t) > for all x and t. 

There are two stationary homogeneous states, i.e. u = and u = 1; the 
latter is stable while the former is unstable against small perturbations. 9 

7.1 Asymptotic behaviour of solutions 

It is well known - and it was proved by Kolmogorov, Petrovskii and Piskunov 
[26] - that if the initial datum is suitably concentrated, e.g. u(x,0) = for 
\x — xq\ > L (i.e. has compact support) as in the cases studied in [31], or more 
generally u(x,0) < Aexp[—x/L], then asymptotically for t — > oo and x — > co 
the solution is of the form u = f(x,t) ~ exp[— (x — vt)/X], with A = \jDje 
and v = \[AeD. This represents a front of width A travelling with speed v; it 
connects the stable state u = 1 and the unstable state u = 0. 

In discussing the FKPP equation, it is convenient to pass to rescaled coor- 
dinates t — et, x — (y/ e/D)x. From now on we will use these coordinates, and 
omit the tildas for ease of notation. In these coordinates, the FKPP equation 
reads 

u t = u xx + u(l - u) . (7.1) 
As for the asymptotic solution given above, this reads now 

u = f (x,t) ~ A exp[-(z-2i)] ; 

note the front has speed v = 2 and width A = 1. 

9 The FKPP equation was introduced in the thirties to model the evolution of the concen- 
tration u(x, t) of a dominant gene in a system undergoing a logistic growth (with time constant 
t = £ _1 ) and such that between two reproductive cycles the individuals move around ran- 
domly with diffusivity D. The same equation also describes autocatalytic reactions; in this 
case u(x,t) represents the concentration of the products of the reaction itself. 



(7.2) 
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It should be stressed that the f(x,t) or fo(x,t) given above provide the 
solution for x — > oo, i.e. in the region 10 of small u; in this region (7.1) is well 
approximated by its linearization around u = 0, i.e. 

u t = u xx +u; (7.3) 

the ansatz u(x,t) = w(z) := w(x — 2t) takes this into the ODE 

w" + 2w' + w = (7.4) 

for w = w(z), with solution 

w(z) = c x e~ z + c 2 ze~ z . (7.5) 

The /o given above, see (7.2), corresponds to c 2 = 0. This can be character- 
ized in terms of symmetry properties, as discussed in the next subsection. 

7.2 Symmetry properties of the linearized FKPP equa- 
tion, and of its asymptotic solutions 

We will now discuss the symmetry properties of (7.3), (7.4) and of the solutions 
(7.5), and characterize the asymptotic solutions in terms of their symmetry 
properties; it suffices to consider translations and scalings. 

It is convenient for our discussion to consider linear combinations of the 
shifts in the (x, t) coordinates, given by X± = d x ^f (l/2)<9 t ; note that z = x — 2t 
is invariant under X-, and that X + = d z . We also write Xq = ud u - Needless 
to say, [Xq, X + ] = 0. 

Lemma 1. The symmetry algebra of the linearized equation (7.3) is generated 
by the scaling X and by the translations X±. The quotient equation (7.4) 
admits only X as scaling symmetry; it also admits the translation symmetry 
generated by X + , while X- has been quotiented out by passing to the z variable. 

Proof. This follows from direct computations. § 

Let us consider the whole set W of solutions described by (7.5); we will 
denote the solution cie~ z + cize~ z as \c\,C2). Note that W = R 2 , and {c\,C2) 
provide coordinates in W. 

Lemma 2. The propagating front solutions are selected among all those in W 
by their symmetry properties; more specifically, they correspond to an invariant 
subspace of W under the action of the group generated by the vector fields X a 
and X + . 

Proof. The transformation e aX ° maps w into e a w without affecting z; the 
transformation e@ x+ maps z — » e 13 z without affecting w. A general element of 

10 Notc that for x — ► and t — * oo we are in the region of u ~ 1; writing u = 1 — /3 we get in 
this region the approximated equation fit = Pxx + 2/3 — 1. With the ansatz /3(x, i) = /3(x — 2t), 
we get the solution /3(z) = (1/2)[1 + e _z (ci cos(z) + C2 sin(z))], with c; arbitrary constants. 
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the group generated by X and X + is written as g(a,(3) :— exp[aX n + /3X + ] 
and acts on W by 

g{a,P) : \c u c 2 ) - |e a+/3 (ci + /3c 2 ), e Q+ ^c 2 ) := |m(ci + /?c 2 ), M c 2 ) . (7.6) 

In the last equality we have defined fi := cxp(a + /3). Thus the action of g(a, (i) 
on W is given, in terms of the (ci, c 2 ) coordinates, by a matrix 

:= (J *) . (7.7, 

The subspace c 2 = is invariant under this action: in other words, the prop- 
agating front solutions (7.2) are indeed selected by a symmetry property as 
claimed. 



7.3 Asymptotic symmetry analysis of the FKPP equation 

It is immediate to see that the only scaling or shift symmetries of the full FKPP 
equation (7.1) are those, with generators X\ — d x and X 2 — d tl corresponding 
to translations in x and t; these reflect the fact that (7.1) is a homogeneous 
equation. 

The situation is different for what concerns asymptotic symmetries, and in 
particular scaling ones, as we now discuss. 

Lemma 3. Let X be a scaling vector field, such that lim^oo exp(AA) extracts 
the behaviour for large \x\ and t. Let Ao be the FKPP equation, and A\ = 
e AA Ao. Then lim,\-»oo Aa = A* is the heat equation u t — u xx . 

Proof. We consider the most general scaling generator, i.e. a vector field in 
the form (3.1), X = axd x + btd t + cud u . We can always set one of the constants 
(a, b, c) equal to unity (provided it is nonzero); this amounts to a redefinition of 
the scaling group parameter. 

Applying the procedure described in previous sections, with of course A := 
ut — u xx — u(l — u) = the FKPP equation, we obtain at once that 

A A = \ c - b [ut - \ b - 2a u xx - \ b u + \ b+c u 2 ] . (7.8) 

We choose c = b < and a = 6/2. As mentioned above, we can set the modulus 
of one of the constants, say b for definiteness, equal to unity; i.e. b = — 1. With 
these choices, we have 

A x = u t - u xx - \~ x u + \- 2 u 2 . (7.9) 

Note we want to extract the behaviour for large x and t; this corresponds to the 
limit A — > oo. The limit A* := limA^oo Aa is the heat equation u t — u xx = 0, 
as claimed. 
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8 Richardson-like logistic reaction-diffusion equa- 
tions 

We now consider reaction-diffusion equations associated to RL diffusion in the 
same way as the FKPP equation is associated to standard diffusion. 
These are written as 

u t = L[u] + h{u) ■ (8.1) 

here L, see (6.2), is the RL linear operator describing passive transport of the 
field u, while h(u) describes its growth. The logistic growth, which we will 
consider here, corresponds to choosing (note an overall constant could be reab- 
sorbed by a rescaling of u) 

h{u) := u(l-u) . (8.2) 
With this choice, and with (6.2), we get the RL logistic (or RLL) equation 

- 2 - a r (2 -a) 



t 



1 — VOL 



+ 



2x 



+ u (1 - u) . 



(8.3) 



(The case v = 1/2 corresponds to standard diffusion, while for v ^ 1/2 we have 
indeed anomalous diffusion). One is usually interested in solutions with initial 
data u(x,0) which are suitably regular and with compact support. 

8.1 Numerical experiments 

A detailed numerical study of RLL equations (8.3) with such initial data was 
conducted in [30, 31] . We summarize the findings of these numerical experiments 
as follows (see [30, 31] for details): 

• (i) asymptotically for large x and t, the solution is described by a trav- 
elling front with varying speed c(t) and width X(t); the form of this front 
for small u is well described by 

; - c(t) ■ t 



u(x, t) ~ A exp 



A(t) 



(8.4) 



• (ii) the (asymptotic) scaling properties of c{t) and \{t) are described by 

c(t) ~ c • t s , X(t) ~ A • t s (8.5) 
where Co and Ao are dimensional constants; 

• (iii) the scaling exponent S is given by 

S := v + (1/a) - 1 . (8.6) 

Thus we rewrite (8.4) in the form 

x- {c t 5 )t~ 



u(x, t) ~ A exp 



x t d 



(8.7) 



Note that for 6 = 0, i.e. for v = 1 — (1/a), the front travels with constant speed 
and width, as for the standard FKPP equation. 
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8.2 Asymptotic scaling invariance 

We want now to describe precisely the invariance properties of the observed 
asymptotic solution (8.7), in particular for what concerns scaling transforma- 
tions. 

Lemma 4. The scaling invariance of the function u(x,t) given by (8.7) is 
described by the generalized scaling group 

x -> ((I s ) x , 
t —* (j, t , 

u — > [exp ((fj, — 1) K t)] u, 
with \i the group parameter. 

Proof. We are interested in vector fields of the form 

X = axd x + btd t + j3(x,t)ud u (8.8) 

which leave (8.7) invariant. It is immediate to note that, in order to leave the 
functional form of (8.7) invariant - i.e. to have at most a rcdifmition of the 
dimensional constant A - we must require a = b5 in (8.8); moreover, we can 
set b — 1 with no loss of generality: this amounts to reparametrization of the 
one-parameter subgroup generated by X. Thus we are reduced to considering 
vector fields of the form 

X = 5xd x + td t + P{x,t)ud u . (8.9) 

The expression of (3 is obtained by the requirement that f(x,t) defined in (8.7) 
is invariant. Applying (1.2) and recalling ip(x,t,u) = (3(x,t) ■ u in the present 
case, this amounts to 

(3(x,t) f(x,t) + Sx(df/dx) - t(df/dt) = 0. (8.10) 

This is immediately solved: we get (3(x,t) = — (c /A )i; for ease of notation, 
we will write K := c /A , so (3{x,t) = —Kt. Inserting this in (8.9), we finally 
obtain that (8.7) is invariant under the scaling-like vector field 

X = 5xd x + td t - Ktud u . (8.11) 

This can be easily checked using (1.2): with u = f(x,t) given by (8.7), and 
writing as usual (8.11) in the form X = £d x + rd t + <fd u , we get Sf := <p = 
^u x - TU t = 0. 

The one-parameter group generated by X acting on the point (x ,t ,u ) is 
given by 

t(s) = e s t , x(s) = e Ss x , u(s) = exp [K(t(s) - t )] u . (8.12') 

This is indeed the generalized scaling group (note that the scaling in u is t- 
dependent) given in the statement. (} 
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8.3 Partial asymptotic scaling symmetry 

Theorem 2. The scaling vector field (8.11) is not a symmetry of the RLL 
equation (8.3). It is an asymptotic symmetry of the same equation. 

Proof. Let us denote (8.3) by Ao and (8.11) by X; applying the second pro- 
longation Y = of (8.7) on A , and restricting to the solution manifold So 
of Ao, which amounts to substituting for w„ according to A itself, we obtain 

A i: = [Y(A )] So = [(l- a )t a& {x/tf- a ]u xx + 

+ [(l/2)(a - l)(a - 2)t aS - 1 (x/t) 1 -*] u x + (8.13) 
- u(l + K - u + Kut) . 

This is not zero, i.e. X is not a symmetry of (8.3). 

Going further on with our procedure, we apply Y on A l7 and restrict the 
obtained expression to So fl Si, with Si the solution manifold of A$; we obtain 

A 2 := [Y(Ai)} SonSi = [2-4Kt + KH*-a(l-Kt)]ui + 

- [(l + K)(a-2)]u = . { ' 

Needless to say, this has the trivial solution u = 0, which is also solution to A 
and Ai , and the nontrivial solution 

_ (2 -«)(! + *) 

U{t > {2-a)-{A-a)Kt + KH* ' (8 ' 15) 

The latter, as easily checked by explicit computation, is in general not a solution 
to A and Ai. Indeed inserting this into A and Ai we have respectively 

Ao = [(a-2)K(l + K)} ( ^ + (2K - 4 + a)t + (2a - 6)\ 

Ai - [(a - 2)K(1 + Kf] ( ^— ^ . (8.16") 

L J \(KH 2 + (a - AK)t + (2 - a)) 2 j 

For K 7^ 0, — 1, both of these expressions are not zero (unless a = 2). However, 
both of these go to zero like (1/i 2 ), for all a and K, in the limit t — * oo. This 
concludes the proof. 11 {} 



8.4 Solution-preserving map associated to scaling 

Having determined that X is an asymptotic (partial) symmetry for our equation 
Ao, we will now apply our general procedure described in section 5 and look for 
asymptotically X-symmetric solutions to Ao. 

11 We could of course also have worked with the symmetry adapted {a, y; w) variables intro- 
duced in the next subsection, i.e. with the expression (8.18) for X and (8.20) for Ao, reaching 
the same result. 
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The first step consists in passing to symmetry adapted coordinates and ex- 
pressing the equation in these. 

The change to symmetry-adapted coordinates and its inverse are in this case 
given by 



a = t , y = x/t 5 
t = a , x = a s y , 



w = ue 



Kt . 

i 

-Ka 



u = we 

In the "new" coordinates, the vector field (8.11) reads simply 

X = od a , 

and the (obviously X-invariant) asymptotic solution (8.7) is 



w = A exp 



JJ_ 

A 



(8.17') 



(8.18) 



(8.19) 



As explained above, once we identify w as the new dependent variable, i.e. 
w = w(a 7 y), the change of variables (8.17') is prolonged to a relation between 
partial derivatives. With standard computations we obtain 



u x 



(w a — 5{y/a)w y — Kw) e Ka 
{1/a 5 ) w y e- Ka 
(l/a s ) 2 w yy 



3.17") 



Inserting (8.17) in the expression (8.3) for Ao, we obtain the expression for 
the latter in the new coordinates; this results to be 



w G = 



2-a 



a* 



w vv + 



y 



+ a 



(!) 



w y + (K+l)w — e 



- K °w 2 



(8.20) 

where we have written x = a(^ — z/ + 1/a) for ease of writing. 

The expression (8.20) holds for the general map (8.17); however we are 
specially interested in the choice 8 = {v — 1 + 1/a), see (8.6). With this, we 
have x = (2 — a), and finally (8.3) reads 



2-a 



W vv + 



2a 



(!) 



l-a 



y 



a- 



Wy + (K + l) w — e 



- K °w 2 



(8.21) 

Note that in the limit a — > oo, the last term disappears (faster than any 
power in a), and (8.21) reduce to a linear equation. 

8.5 Scaling-invariant (asymptotic) solutions 

Theorem 3. The equation (8.21), has no nontrivial X -invariant solutions. It 
admits nontrivial asymptotically X -invariant solutions. 
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Proof. The X-invariant solutions to (8.21) are obtained by requiring that 
w a = 0; with this the equation reduces to 



\a 



2-a 



W v 



2-a 
2a 



(!) 



l-Q 



u 



w y + (K + l)w - e- Ka w 2 



. 



(8.22) 

Note that a appears parametrically here, and (8.22) splits into the equations 
corresponding to the vanishing of coefficients of different powers of a (this is a 
general feature of partial or "weak" symmetries, see [10, 39]). The only common 
solution to these is w = 0, which ends the proof of the first part or the statement. 

We go back to considering (8.21); in order to study its asymptotic behaviour 
for a — > co, we disregard the term which is exponentially small for large a. The 
resulting linear equation for w — w(y), i.e. 



2-a 



2-a 
2a 



(!)' 



y 

a— 
a 



w y + (K + l)w = , (8.23) 



(8.24) 



(!) 

yields as solution 

w(y) = d/C[0,0] + ^2(K + l)a^y/ar) C 2 K[l/2, 1] 
where C\,C<2, are arbitrary constants, and 

K\x,y\ := F n [(K + l)a/a 2 + x; 1/2 + y; -a(y/a) a 
with Fn = \F\ the Kummer confluent hypergeometric function 

r(6) 



Fu[a;b;z] := iFi(a;&;z) 



r(6-o)r(a) Jo 



e zt t a-l 



(1-t) 



b-o-l 



(8.25) 



At. (8.26) 



Needless to say, the asymptotic solution (8.24) could be expressed in terms 
of the original variables using (8.17); this yields an explicit but involved and not 
specially illuminating expression. ^> 



9 Other asymptotic partial symmetries of RLL 
equations 

In the previous section we have analyzed the RLL equation (8.3) on the basis 
of the scaling symmetry (8.11) of its asymptotic solution (8.7). However, this is 
not the only symmetry of the observed asymptotic solution (8.7). In this section 
we apply our approach on the basis of different symmetries. 

9.1 General symmetries of the asymptotic solution 

Lemma 5. The vector field X = £d x +Tdt+<pd u is a symmetry of the asymptotic 
solution (8.7) to the RLL equation (8.3) if and only if it belongs to the two 
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dimensional module (over smooth real functions C°° (R 3 , R) ofx,t,u) generated 
by 

Xl = dx+ ( a^f) 9 « J (9 i) 

X 2 = (X t 1+S ) 8 t - (xS + c t 1+s ) 8 U . 

Proof. This follows easily by using (1.2) and the explicit expression (8.7) of 
the asymptotic solution u = f*(x,t). Indeed, applying (1.2) we get 

^ ~ ^ = ip - A exp [(cot - xt- s )/X ] (x5 + c a t 1+s )r - t£ overX t 1+s , 

and the result follows immediately ()■ 
In the following, we will consider in particular 12 

X := (xS + c t 1+s ) 8 X + td t , (9.2) 

as well as X\ and X 2 themselves. 

We will write second-prolonged vector fields in the form 

Y = X^ = X + * x A + * t ^_ + $ xx ^_ + * It ^_ + * tt ^_ ; (9.3) 

ou x du t du xx du xt ouu 

we will actually need only the coefficients \I> X , \I> t , ^f xx . 

We will, as in the previous discussion, denote (8.3) as A . Let us start by 
considering X±. 

Theorem 4. The vector fields Xo, X\ and X 2 are partial symmetries of Aq. 

Proof. We will denote by Yj the second prolongations of Xi. Let us start by 
considering X x ■ In this case the coefficients of the second-prolonged vector field 
Yi are: 

* x = , tt t = -(d/X Q )t-^) ,* xx = 0, * xy = 0, *„ = (5/X )(l+6)t- (2+ V . 

(9.4) 

We define then 

Ai = [Fi(A )] So ; 
A 2 = [Y^A,)}^ . 

By explicit computation, it results that 

[^i( A 2)] 5onSinS2 = ; 

this shows that Xi is a partial symmetry for A . 

Let us now consider X 2 . The relevant coefficients of Y 2 are 

8 t (1 + 8)5x „, 



(9.5) 



12 Note that u is not acted upon by this, i.e. we have a vector field which is horizontal for 
the fibration (M,ir,B). 
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We get more involved expressions for Ai and A2 obtained as above (we mention 
that in this case, A2 identifies a single function u = u(x, t)), but again 

[F 2 (A 2 )] SonSinS2 = ; 

this shows that X 2 is also a partial symmetry for Ao. 

Finally, we consider the vector field Xq, in which we are specially interested. 
The relevant coefficients of Y are 

= Su x , * t = -(1 + 5)c t s u x - u t , * xx = -2Su xx . (9.7) 

Defining Ai — Y (A ) and A 2 = lo(Ai), we obtain again that 

[^o(A 2 )] 5onSinS2 = , 

hence Xq is also a partial symmetry for Ao. ^> 

9.2 Invariant solutions 

It turns out that reduction, and invariant solutions, under the vector field X 
are of special interest. This is due to the following theorem, which provides an 
analytic explanation of the numerically observed behaviour. 

Theorem 5. The RLL equation (8.3) admits an asymptotically X^-invariant 
solution, described by (8.7). 

Proof. In this case the symmetry adapted coordinates are 

a = t , y= (x/t 5 ) -c t , w = u . (9.8) 
The inverse change of coordinates is therefore given by 

t = a , x = (y + c o-)a s , u = w (9.9) 
and the relevant u derivatives are expressed in the new coordinates as 
f Sy + cp(l + S)ct \ 1 1 

U t = W a - [ ) Wy , U x = JgWy , U xx = -pgWyy . (9.10) 



u 



The equation (8.3) is therefore written, in these coordinates, as 

Wa = ( coO+^W j Wy + [ (2 - a) Wy + 2(c o- + y)w vv ] 

+ (1 — w)w . 

(9.11) 

In these coordinates X reads simply 

= od a , (9.12) 

and its second prolongation is 

d d d d 

Yq = a- w a - w ay - 2w TO - . (9.13) 

oa ow a ow ay ow aa 
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In order to discuss (9.11), it is convenient to rewrite it as 

W a = AlVyy + BlVy + f (w) , 

with 



(9.14) 



a 



2-a 



B 



c + o h e 



a 



V + CqCJ 
a 



2-a 



1 



For g — > co , (9.14) reads 



c a w yy + c {l + S)w y + f(w) . 



y + c Q a J 
(9.15) 

(9.16) 



The Xp-invariant solutions satisfy w a = 0, i.e. w = w(y), and are thus obtained 
as solution to 

^"Wyy + C (l+6)Wy + W = , (9.17) 

where we have used w « 1 in the region we are investigating, i.e. for a — > oo, 
so that f(w) ~ w. 

Solutions to (9.17) are of the form 



w ( y ) = Cl e-"+ y + c 2 e-"-y , 



(9.18) 



where u>± — ^t S i 



l±Jl 



. If we require the solutions to be non 



oscillating, this implies a lower bound on the parameter Co, i.e. cq > (2/(1 + 
5)f' a . 

The solution e~" +z is unstable against small perturbations, while e~ u - z 
is stable [26]. As proved by Kolmogorov [26], the asymptotic solution is the 
stable one with the lowest speed giving nonoscillating behaviour, i.e. cq = 
[2/(1 + 8)] 2 ' a . This means w(y) ~ with u = [2/(1 + S)} 1 - 2 ^. Going 

back to the original variables, we get 



u{x, t) ~ A exp 



x — v(t)t 



A exp 



-ujq- 



c t 1+s 



(9.19) 



This is precisely the numerically observed asymptotic behaviour, described by 
(8.7); see [31]. ' 



10 Asymptotic symmetry in an optical lattice 

In this section we apply our approach to a different kind of anomalous diffusion 
equations; that is, we focus on the equation describing anomalous transport 
in an optical lattice [29]. The equation for the marginal Wigner distribution 
w(p, t) of the momentum p at time t reads 

d 

w t = -j- [h(p)w - g(p)w p ] (10.1) 
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where the functions h(p) and g(p) are given by 

m := T+fe 5 ' 9(P) := 70 + TT(kF ; (10 ' 2) 

here a, 70, 71 are certain constants, and p c represents the capture momentum. 

With the shorthand notation j3(p) := 1/[1 + (p/p c ) 2 }, the equation (10.1), 
(10.2) reads 



w t = (7o+/3(p)7i) w pp + (af3(p) - 2-/ 1 [(3(p)/p c } 2 ) pw p + 
+ a (/?(p) - 2[(p/p c )/3(p)] 2 ) W . 



(10.3) 



Attention to eq. (10.3) was recently called by Lutz [29] (to which the reader 
is also referred for derivation and a discussion of this equation), who remarked 
that - quite surprisingly the equilibrium distribution is related to Tsallis 
statistics [43]. Indeed, for q < 3 the stationary solution of (10.3) is the Tsallis 
distribution 

Mp) = \ [1 - (/?//*) p 2 ]" (10-4) 

where 



P = — , I x . A* = q ; , 9 = 1 + := 1 + <5 



— =1 + ^ 

2(70 +7l) 1 ~ 1 a Pl 

and Z is a normalization factor, which for g < 3 can be chosen so that wo(p)dp 
1. For 5/3 < g < 3 the second moment / p 2 'u;o(p)dp of the Tsallis distribution 
is infinite, and we have anomalous diffusion. 

Here we will consider generalized scaling transformations 13 ; that is, trans- 
formations acting as a standard scaling in the independent p and t variables, 
and as a p and t dependent scaling in the dependent variable w(p, t): 

X = cipdp + c 2 td t + tp(p,t,w)d w . (10.5) 

We will of course consider asymptotic invariance. 

Applying (1.2) to (10.4) and to generalized scalings we get with easy com- 
putations that (up to a common factor) the transformation generated by (10.5) 
with ci = — 1 leaves (10.4) invariant if and only if 

if = (2f3Z q - 1 )p 2 w q ; (10.6) 

obviously C2 remains unrestricted at this stage. We will, for ease of notation, 
write 

v := 2[3Z q - 1 . 
Let us now discuss what is special in the vector fields 

X = -pd p + atd t + [vp 2 w q ]d w (10.7) 



13 The peculiar properties of (10.3) and its stationary solution, in particular concerning the 
scale-free nature of wo(p), has been studied by Abe [1] in connection with dilation symmetries 
and canonical formalism. Note however that he introduces a canonical structure, which wc 
do not need here, and that his symmetry generators are nonlocal, being expressed as integrals 
over p. Thus, our approach is substantially different from his one. 
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for what concerns their action on the equation (10.3). 
We are interested in 

Y = X + * t A + + * pp 



dw t dw p dwpp 



(other prolongation coefficients are irrelevant for application to our equation); 
the coefficients corresponding to (10.7) are 



2 w s ; 



\fr t = -ow t + vqp 2 w s ; * p = w p + 2upw ( - 1+s '> + v(l + 5)p 
* pp = 2w pp + 2vw {1+s ^ +vw( s -^ [6(1 + 6)p 2 w 2 + (1 + 5)pw(Aw p + pw pp )] . 

Applying Y on (10.3), and substituting for w t according to (10.3) itself, we 
obtain an expression which we write as 

A = T (p,t,w) + T 1 (p,t,w)w p + T 2 (p,t,w)w pp . (10.8) 

The explicit expression of T 2 is 

P 70 (2 + a) (Pc 2 + p 2 f + 7i Pc 2 (Pc 2 (2 + a) + (4 + a) p 2 ) 

1 2 = 2 ' 

(p c 2 +p 2 ) 

The limit of T 2 for large \p\ is nonzero unless we choose 

a = -2 . (10.9) 
We assume this from now on. With (10.9), X reads 

X = -pd p - 2td t + vp 2 w q d q (10.10) 

and T 2 reduces to 

-2 llP 2 p 2 

A 2 



(Pc 2 +P 2 ) 2 



therefore, r 2 — > for \p\ — > oo. 
As for Ti, with (10.9) it reads 

ri = [2p((p^+P 2 )(a^-2(l + J) 7 o^(^+^) 2 ) 
-27ip2((l + 6)vp^w s +p 2 (-l + vw d p 2 + 6vw 5 p 2 ) 

+p 2 c (l + 2vw 5 p 2 + 25vw 5 p 2 )] x [(p 2 c +p 2 ) 3 ] -1 . 

It can be checked that in the limit \p\ — > oo and to — > 0, we get Ti — > 0. 
Finally, with (10.9) T reads 

To = [ W (-2^ 5 (p2 + p2 ) ( 7lp 2 (:p 2_ p 2 )+7o(:p 2 +:p 2 ) 2 ) + 

+a;> 2 (-((2 + 6)uw s p 6 ) - 2p2p2( 3 + 2vw s p 2 )+ 
+p i c (2-2vw 5 p 2 +6vw s p 2 )))] x [p 2 c + P 2 y 3 . 

Once again, in the limit |p| — > oo and iD^Owe get T — > 0. 
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This explicit computation suggests that 

X = -pd p -2td t + [(20Z"- 1 )p 2 w' 1 ]d w (10.11) 

is an asymptotic symmetry for the equation (10.3), provided this is set in the 
function space identified by \imi p i ^ ao w(p,t) = 0. 14 

This limitation is removed by a more careful analysis, using the systematic 
procedure described above; we are now going to describe this approach. 

Symmetry adapted variables for X are given by 

y = p 2 /t , v = w~ 5 - {v5/2)p 2 , a = t; (10.12) 

the inverse change of coordinates is 

t = a , p=jya , w = (v + {u8/2)y<j)- 1 ' s . (10.13) 

In the new coordinates the vector field (10.11) is simply X = —2crd a . Its second 
prolongation is Y = X + 2v a {d/dv a ). 

With standard computations (see sect.l), the partial derivatives appearing 
in (10.3) are written in the new coordinates as 

W t = W a ~ {y/v)Wy 

= -{l/{5a))[{5vay/2)+v]- q/5 [av a - yv y ] , 
w p = 2^/y/aWy 

= -{2^/8)[(5 V ay/2)+v]- q/S [(5u*/2)+v y ] , 
w P p = {2/a)w y + (Ay/a)wyy 

= 2 [-25v(Sua + 2v y + Ayv yy )+ 

+y(25(4 + ZS)vov v + 8(1 + 5)v 2 y + 5 2 va{(2 + S)va - 4yv m )j\ x 
x [5 2 a{8vay + 2v) 2 ]^ . 

(10.14) 

Using these we can rewrite the equation (10.3) in X-adapted coordinates; 
this is a rather involved expression, and we write it as 

A = (xiv<7 + X2V yy + X3^l + XiVy + 9) (v+ {\/2)5vayY 1/S , 

14 More precisely, we should also check that the variation Sw = p 2 w q (p,t) is normalizable. 
This is a more restrictive requirement, but still verified for functions near to the Tsallis 
distribution (10.4) for all q in the physical range 1 < q < 3. However, this limitation will be 
removed below. 
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where Xi an d are function of a, y, v alone, given by 
Xi = -2 [d(2v + Suay)}' 1 , 

X2 = 8y ((70 + 7i)^ + 7 cry) [&7-(#! + <tj/)(2v + Suay)] 
X3 = l6y{l + 5)((-fo + "f 1 )p 2 c +"f (7y) [S 2 s(p 2 c + (jy)(2v + Siyay) 2 y , 
Xi = -2 (-y(p 2 c + °y)(p 2 c - 2ap 2 c a + ay)(2v + 5uay)+ 
+ 2 l0 ( P 2 c + <7y) 2 (-2v+(4 + 3S)v<jy)+ 

+ 2^ip 2 (-2p 2 v + Avp 2 ay + 35vp 2 c ay + 2avy + Ava 2 y 2 + h5va 2 y 2 )) x 
x Sa(p 2 + ay) 2 (2v + 5vay) 2 \ 1 
9= -[2^v(p 2 c +ay) 2 {-2v + (2 + 5)vay) + 

+ p 2 (a{2v + 5vay) {—2p 2 v + 2vp 2 c ay — b~vp 2 c oy + 2avy + 2vo 2 y 2 + 8va 2 
+ 2^]v{-2p\v + 2vp 2 c oy + 5vp 2 c ay + 2avy + 2va 2 y 2 + 35va 2 y 2 ))] x 
x [(p 2 c +ay) 2 (2v + 5u<jy) 2 ]- 1 . 

Applying Y on (10.3), we obtain a quite complex expression, which we write 



as 

Ai = y(A ) = [E 1 v yy + E 2 v 2 y + E 3 v y + e] (v + (1/2) 5vay)- 1/S (10.15) 

where Sj and do not depend on v y , v yy . The explicit expressions of these 
are not simple, so we will provide only their scries expansion around s = 00 in 
terms of e — 1/s. These are: 

51 = 16 [((1 + 2S)g ) I (vP)\ e 2 + 0(e 3 ) ; 

5 2 = - 32 [(1 + 5)(1 + 3,5)30 / (v 2 S 5 y)] e 3 + 0(e 4 ) ; 

~ 3 = 4 [((1 + 2S)[S(y - 2ap 2 c ) - 2(4 + 35)}g ) / (vS^y)] e 2 + 0(e 3 ) ; 

6 = 2 [((2 + 35 + 5 2 )(2g + a5p 2 c )v + 2(1 + 25) v a ) / {v5 3 y)} e + 0(e 2 ) . 

(10.16) 

Needless to say, these all vanish for e — > 0, i.e. for s — > 00. This computations 
shows that: 

Theorem 6. The one-parameter group of transformations generated by 

X = -pd p - 2td t + [(2l3Z q - 1 )p 2 w q ]d w 

is not an exact symmetry of the equation (10.3), but is an asymptotic symmetry 
for this equation. No vector field of the form (10.7) with a ^ — 2 is an asymptotic 
symmetry of (10.3). 

In other words, the distribution (10.4) is an invariant function for a trans- 
formation X , given by (10.10), which is an asymptotic symmetry of (10.3), and 
hence the scale-free nature of the asymptotic solution (10.4) is a consequence of 
the asymptotic symmetry properties of the equation. 

11 Discussion and final remarks 

Let us briefly summarize the main findings of previous sections. 
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In sections 1-2 we recalled some basic facts about geometry of PDEs and 
the concepts of full, conditional and partial symmetries of these; in sections 
3-4 we developed our approach, based on geometrical considerations, to asymp- 
totic versions of these. In sect. 5 we described a method to test the asymptotic 
behaviour of nonlinear PDEs based on asymptotic symmetry properties. The 
rest of the paper is devoted to application of the general method described and 
developed in sects. 1-5 to some specific classes of PDEs. 

First of all we considered, in sect. 6, a model class of anomalous diffusion 
equations, which we called "Richardson-like" (RL), see (6.2), which had been 
previously studied numerically - in particular for what concerns their asymptotic 
behaviour - in detail [31]. With theorem 1, we proved that the asymptotic 
properties can be recovered by the standard ones for the heat equation, via the 
change of variables (6.3), and are described by (6.6). 

In sect. 7 we applied our approach to the standard FKPP equation, and 
described in detail its asymptotic symmetry properties, see lemmas 1-3; we have 
also shown that our approach recovers the well known asymptotic properties of 
FKPP solutions. 

In the following two sections we have then applied our method to anomalous 
reaction-diffusion equations associated to RL type, also called "Richardson-like 
logistic" (RLL) equations, which were also studied numerically in [31]. In sect. 8 
we recalled the features of asymptotic solutions as observed in numerical ex- 
periments, and identified the Lie generator X, given by (8.11), of the observed 
asymptotic generalized scaling invariance (lemma 4); in theorem 2 we showed 
that this is not a symmetry of the RLL equation, but it is an asymptotic sym- 
metry for it. We have then considered the solution-preserving maps associated 
to this asymptotic scaling symmetry, focusing on the physical value of the pa- 
rameter 5; in theorem 3 we have shown that in this case our equation has no 
solution invariant under X (which therefore is not a conditional symmetry of 
the equation), but admits solutions which are asymptotically invariant under it: 
X is an asymptotic conditional symmetry of the model anomalous RD equation 
we consider. 

In sect. 9 we passed to consider in more detail the numerically observed 
asymptotic solution (8.7). With lemma 5 we identified the full symmetry alge- 
bra Q of it; vector fields in this algebra are asymptotic conditional symmetries 
of our model RD equation. In theorem 4 we focused on certain vector fields in 
Q, and shown they are partial symmetries for the model RD equation. Among 
these vector field is the scaling vector field Xq given by (9.2); it does not depend 
on, nor acts upon, the dependent variable u. With theorem 5, we proved that 
the RLL equation does admit an asymptotically X -invariant solution, which is 
precisely the observed one (8.7). 

Finally, in sect. 10 we consider a different example of anomalous diffusion, oc- 
curring in optical lattices. This attracted some attention due to the appearance 
of a Tsallis distribution as stationary states in a certain range of parameters. 
The surprising fact is that such a distribution decays polynomially, and there- 
fore is scale-free. Our method can deal successfully with this equation, and 
indeed theorem 6 identifies the relevant asymptotic symmetries. Moreover, the 
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observed (scale-free) properties of the solution in the physically relevant range 
can be explained on the basis of the asymptotic symmetry properties of the 
underlying equation. 

Final remarks 

Let us present some final remark. 

1. ) One could wonder what is the advantage of the discussion conducted 
in sect. 9, and in particular of theorem 5, over the one of sect. 8 and theorem 
3. The answer is in the properties of Xo: it identifies characteristics (in the 
sense of the elementary theory of linear or quasilinear PDEs) in the space of 
independent variables (x,t), and reconducts the observed asymptotic properties 
of solutions to the RLL equations to their asymptotic constancy along these 
characteristics. This feature should not be surprising: indeed we had already 
observed in sect. 8, cf. eqs. (8.22) through (8.23), that the asymptotic regime is 
described by a quasilinear equation. Our result could then be also described in 
terms of " asymptotic characteristics" . 

2. ) Our general method to investigate and predict asymptotic symmetry - 
in particular, scaling and front-like - properties of solutions to scalar PDEs, 
represents an evolution of the classical method to determine partially invariant 
solutions for symmetric PDEs [37, 45], and a blend of it with the method of 
conditional and partial symmetries [12, 24, 28, 45], in order to analyze equations 
which do not have complete (as opposed to asymptotic) symmetries - as indeed 
for the equations considered in this note. Our method is based on the abstract 
approach developed in [20], based itself on ideas and previous work by several 
authors [4, 7, 14, 23]. 

3. ) The main body of this work was concerned - and the method described 
here can strictly speaking deal only - with scale invariance (at infinity or near 
a travelling front); this is surely not the most general kind of (asymptotic or 
exact) invariance, but definitely one of physical relevance. We trust that suitable 
generalizations of our approach can also deal with more general asymptotic 
invariance properties. 

4. ) Here we focused on a restricted class of anomalous diffusion and reaction- 
diffusion equations, which one of us had previously studied numerically [31], and 
on the equation for the marginal Wigncr distribution in an optical lattice, whose 
stationary solutions are known to have certain asymptotic scaling properties 
[1, 29]. It is however quite clear (also due to the geometrical nature of the ideas 
at the basis of our constructions) that the method developed and applied in 
this paper does apply to quite general differential equations, and is potentially 
capable to provide a sound explanation - or prediction - of the asymptotic 
invariance of their solutions. 
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